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Abstract 

We determine the holographic spectra of scalar mesons from the fluctuations of the em- 
bedding of flavor D-brane probes in HQCD models. The models we consider include 
a generalization of the Sakai Sugimoto model at zero temperature and at the "high- 
temperature intermediate phase", where the system is in a deconfining phase while ad- 
mitting chiral symmetry breaking and a non-critical 6d model at zero temperature. All 
these models are based on backgrounds associated with near extremal Nc D4 branes and 
a set of Nf « N^. flavor probe branes that admit geometrical chiral symmetry breaking. 
We point out that the spectra of these models include a branch which does not 
show up in nature. At zero temperature we found that the masses of the mesons Mn 
depend on the "constituent quark mass" parameter mj^ and on the excitation number n 
as ~ and ~ n^''^ for the ten dimensional case and as M„ ~ and M„ ~ 7^0.75 
in the non-critical case. At the high temperature intermediate phase we detect a decrease 
of the masses of low spin mesons as a function of the temperature similar to holographic 
vector mesons and to lattice calculations. 
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Whereas realizing confinement in dual holographic models of QCD (HQCD) is easy 
the incorporation of flavored chiral quarks and in particular chiral symmetry breaking is 
more difficult. Sakai and Sugimoto [1] proposed a model that admits the two phenomena. 
It is based on placing a set of Nj D8 and anti D8 probe flavor branes into the gravity 
model of near extremal D4 branes [21 13] . 

The mesonic spectra is one of the most important properties of hadron dynamics that 
can be "measured" in the HQCD laboratory. The low spin mesons are associate with the 
fluctuations of the fields that reside on the probe flavor branes, the vector mesons with 
the U{Nf) flavor gauge fields and the scalar mesons with the embedding of the probe 
branes 0. Here in this paper we focus only on scalar mesons. The motivation behind 
addressing this problem are the following: (i) To verify that the meson spectrum at zero 
temperature does not include tachyonic modes. Had there been such modes it would 
have indicated that the system is unstable. Since the model of [1] is based on placing 
branes and anti-branes one may be worried that the system is unstable and hence the 
importance of this check, (ii) The spectrum of the scalar mesons has been determined 
already in [1]. However the attemps to derive it in generalizations of the model where the 
asymptotic separation of the brane anti-brane L is smaller than half of the circumference 
of the compactified direction X4, namely for L < irR failed for the symmetric modes [5],[6] 
(iii) To determine the dependence of the spectrum on the excitation number n and the 
parameter defined in ( H2|) that is related to the constituent quark mass. In addition 
one naturally would like to compare the explicit ratios of meson masses that one deduces 
from any given HQCD model and the experimental data to get an indication of how 
well the model describes real hadron physics, (iv) To further examine the differences 
of physical properties extracted from critical models to non-critical models which were 
previously discussed in O [3, [8] . The spectrum of scalar mesons was extracted also in 
other HQCD models[9l [TOl [11]. For further reading see [12] and references therein. 

We can summarize the outcome of the paper as follows 

• We were able to choose coordinates that avoid the singularities that were en- 
countered in previous works [3 [6] and determine the spectrum of both the anti- 
symmetric as well as symmetric branches. 

• We find that in the models examined and in particular the original model of [1] the 
symmetric solutions correspond to scalar mesons of the form 0"'"'' whereas the anti- 
symmetric solutions correspond to mesons. This property which seems to be 
in common to a HQCD models based on probe branes and anti-branes, contradict 
the low lying spectrum in nature. There are no low lying mesons. 

^ High spin mesons are naturally described by semi-classical spinning string configurations [3| 



1 



• At zero temperature we found that the masses of the mesons Mm depend on the 
"constituent quark mass" and on the excitation number n as M^^ ~ and 

~ with a ~ 1.7 for the ten dimensional case and as ~ and ~ 
with (3 ~ 0.75 when a CS term is incorporated and ~ 1 without such a term 
in the non-critical case. At the high temperature intermediate phase we detect a 
decrease of the masses of low spin mesons as a function of the temperature similar 
to holographic vector mesons and to lattice calculations. 

The paper is organized as follows. We begin in section [T] with a brief review of the 
holographic models we investigate. We summarize the main features of the model of 
Sakai and Sugimoto at zero and finite temperature and an analogous six dimensional 
non-critical model. In section [2] we describe the extraction of scalar mesons from the 
fluctuations of the embedding. In particular we point out that in the coordinates in- 
troduced in [1] the eigenvalue problem admits a singularity that prevents the numerical 
determination of the eigenvalues. A different coordinate system is presented in section 
|3] which evades the problem of the singularity. Using these coordinates, the spectrum 
of masses as a function of the constituent mass and excitation number is derived. The 
spectrum of scalar mesons that follows from a non critical model of near extremal D4 
branes is analyzed in section HI Section [5] addresses the issue of parity and charge con- 
jugation of the scalar mesons. It is pointed out that the spectrum includes mesons 
which do not show up in nature. Section [6] is devoted to the spectrum of mesons above 
the deconfining phase transition in the "intermediate phase" . We summarize the results 
and raise certain open questions. 



1 Review of the holographic models 



1.1 The Sakai Sugimoto model 

The model of [2j, describes the near horizon limit of Nc Z)4-branes wrapping a circle 
in the X4 direction with anti periodic boundary condition for the fermions. Into this 
background a stack of Nf D8 is placed at ^4 = and a stack of Nf D8 at the anti-podal 
point of the circle [1]. Assuming Nf << Nc one can overlook the modification of the 
metric and dilaton due to the backreaction of the background by the Nf D8-D8 systems 
and continue to use the metric and dilaton associated with the Nc D4 alone. Therefore 
the metric, dilaton and the RR four form are given by 
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Where V4 denotes the volume of the unit sphere and €4 its corresponding volume form. 
Ig is the string length and Qg a parameter related to the string coupling. The X4 is the 
compactified direction that is asymptotically transverse to the D8. The manifold spanned 
by the coordinate u, has the topology of a cigar where its tip is at the minimum value 
of u which is u = ua- The periodicity of this cycle is uniquely determine to be 

3:,^='^(flm\ =2nR (2) 



3 V ua 

in order to avoid a conical singularity at the tip of the cigar. The classical profile of the 
D8 probe brane in this background is given by the classical solution to the e.o.m of the 
DBI action of that probe brane. The D8 DBI action is 



Sns =TsJ dtd'xdudn,e-^4^g = f,j dtd^xduu^ 
= Ts J dtd^xdx^u 



+ ^ (3) 



«ov/(%) = corist (5) 



where g stands for the pullback metric on the D8 brane. The simplest way of solving 
this e.o.m is by noting that the action is independent of X4 and so its Hamiltonian is 
conserved. 

u'fju) 

^/(■") + (^)'7& 

where we assumed that there is a point uq where the curve u{x4), which describes the 
profile of the D8 brane in the (m, X4) plane has a minimum. At that point the D8 brane 
bends, namely the D8-D8 join together. After some algebra one finds 

fdx4\ 1 

(6) 



Hence we find that the profile of the D8 brane probe is 

Where uq is a constant of integration setting the lowest value of u to which the D8 
barne is extending. At that point the D8 brane join the D8 brane and the brane is 
extending back into the UV. The value of Uq also sets the asymptotic distance L between 
the position of the D8 and D8 brane 

dx4 = 2r- = 2(^) r dy 'T^ (8) 
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Hence we see 



' \ 1/2 

Loc(^) (9) 
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For later use we define 



u 

^ = f{u)u^ - fiuoK 

The DBI action then becomes 

S = TsJ e-*^J\detgo\ ~ J d^xdu-i^'^u^'^ 

1.2 Thermodynamics of the Sakai Sugimito model 

In [I3]|!| a study of the thermodynamics of the Sakai Sugimito model was carried using 
the conjecture presented in [2]. 

The conjecture states that the thermodynamics of a field theory with a gravitational dual 
is determined by taking into account the contribution to the saddle point approximation 
from all the gravitational backgrounds with the correct 'UV asymptotic, with compacti- 
fied Euclidean time direction of period (3 = ^ and with anti-periodic boundary condition 
for the fermions along this direction. The temperature of the field theory is T = and 
its properties are read from the manifold responsible for the most dominant contribution 
to the saddle point approximation namely the one that has the lowest free energy. 
When ever one background looses its domination to another background as we vary the 
temperature, a phase transition occurs in the dual field theory. 

In [12] two manifolds where found to have the same 'UV asymptotic as the one of 
Sakai and Sugimoto model, the background ([T]), and the same configuration only with 
the time and X4 directions interchange. 



ds^=[^][-f{u)de+ Si.dx'dx^ + dxl] + (— 
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/(«) ^ 1 - (^)' (12) 

and the temperature is given by 



St = ^r-^ =(3 (13) 
3 V Mr 



^see also [H]. 
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The difference between the free energy densities of the two backgrounds is proportional 

This means that when the circumference of the X4 cycle is 
smaller than that of the time direction namely when T < 1/27: R the background ([1]) 
is the dominant one, while when the opposite occurs and T > 27iR the action of (ITTil 
will dominates. At the temperature T = 1/2ttR the two actions are the same since the 
two backgrounds are different by the labeling of the coordinates, so at T = Tc = 1/2'kR 
the system has a first order phase transition. In [13] it was argue that in the dual field 
theory, the physical interpretation to this phase transition is a transition from confined 
phase at T < 1/2'kR to deconfined at T > 1/2ttR. This can be seen via a computation 
of the quark anti-quark potential [15] in the two backgrounds. Another indication to this 
interpretation is that the renormalized free energy of the low temperature phase shows a 
A^^ behavior while that of the high temperature phase shows a one. Hence from now 
on we will denote Tc = T^. 

At the high temperature phase there is another possible classical solution to the profile 
of the D8 brane which is a configuration with constant X4 namely, 2:4 (m) = 0, 
Now since the bulk free energy is the same for the two configurations of the D8 branes, the 
difference of the free energy of the D8 probes determines which of the two configurations 
is the preferable one for a given temperature. It turns out that the transition between the 
two configuration depends on the parameter ?/t = its value at the phase transition 
turns out to be ~ 0.73572. 

Using eq.® we find Lc = 0.751 ( j , hence at the critical point i/t = Vt the critical 
temperature is set by the asymptotic distance between the branes (setting Rd^ = 1) 



Tc = -^4^' = ^{y^T^oy^' = 0.154/L (14) 

The field theory sees this transition as chiral symmetry restoration at high temperature, 
this interpretation is natural since the D8 branes are now disconnected and there is an 
U{Nf) X U{Nf) global symmetry. 

Hence we will denote this critical temperature as T^sb- Note that this only happens at 
the high temperature phase so there is still the condition T^sb = 0.154/L > l/27iR. 
So if L > 0.97-R, we find that Ta is always higher than T^sb, and so deconfinement and 
chiral symmetry restoration phase transition happen together. We see that in this model 
xSB and confinement appear independently of one another as a result of the existence 



Of course in our model there are also D8 brane which their DBI action will contribute to the 
total free energy of the configuration as well, but this is sub-leading to the bulk action since the bulk 
action is of order and the contribution of the D8 is of order Nc ■ Nf which is negligible in the probe 
approximat ion 

^This configuration was not possible in the low temperature, but in the high temperature phase the 
time circle shrink to zero at u = ua and so the D8 brane can just smoothly end there. 
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of the free parameter L coming from the hd nature of the field theory. 
1.3 Non critical holographic model 

A non critical model with a very similar properties to Sakai-Sugimoto model was pre- 
sented in [5l [7] 0, this model consists of non-extremal configuration of Nc DA branes 
placed in a six dimension space-time with one of the DA coordinates taken to be periodic 
with anti periodic boundary condition for the fermions. 

The metric, dilaton and RR six-form field take the form [7] 




with 

/H = l-(^)' (16) 

and where Qc is proportional to Nc-, the number of color DA branes. In order to avoid 
conical singularity the periodicity of the cycle of X4 is set to 

0:4 ~ 0:4 (5x4 ; (5a;4 = ^'^^^'^^ (17) 

Of course the curvature of order one of this background makes the leading order super- 
gravity an un justified approximation to string theory on this background. Nevertheless 
its believed that at least the extremal model due to its symmetries, is indeed a good 
background for the study of non-critical string theory [16]. Now we place Nf DA branes 
which are transverse to the cycle and extend up to infinity in the u direction. The 
properties of the four dimensional low energy effective field theory living on the intersec- 
tion of these color and flavor DA is then seem to be very similar to those found at the 
Sakai Sugimoto model. Thus we would like to study its spectrum of scalar excitations 
and check if there is no tachyon in the model. 

Just like in the critical model the DA brane may bend on the (m, X4) cigar and in order 
to find its profile one must solve the e.o.m of the X4 coordinate. This e.o.m is derived 

^For other non-critical SUGRA models with flavor see [TCI [T71IT51 [1^1 . 
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from the action of the flavor DA brans namely the DBI action plus the CS term which 
are are given by 

= -T4 J d''xe-'^^-det{g) + T^a j P{C^^^) (18) 
Following similar steps to does taken in the previous section we find 

x,,iu) = r iuimuo)-aul + au-)du' ^^^^ 

where 0, = 

2 Fluctuation of the embedding and scalar mesons. 

We now turn our attention to the study of the fluctuation of the D8 brane around its 
classical profile. As was mentioned in the itroduction, one has a two fold interest in 
those fluctuations: (i) They correspond to scalar mesons in the dual gauge theory, (ii) 
Tachyonic modes of the fluctuation signals an instability of the system. 

We start by expanding the coordinate around its classical value and deflne the 
fluctuation ^{u,x^) as follows: 

X4(m, x'^) = X4(M)d + i{u, x^") (20) 

Substituting this into the action ([3]) and expanding to quadratic order in ^ we flnd the 
following action for the fluctuations 

So,\j d'xdu[u"'R],,r"'v"'d,id.i + u'"'r"\d^if] (21) 

where 7 is deflned in ( ITUl) . We now introduce the following mode expansion 

00 

i{u,x^) = Y.Ux^)Uu) (22) 



n=0 



Using the symmetries along the x^ directions we have 

V^'d^dJ^ = -mlU (23) 
The e.o.m for the ^„ modes reads 

d4{u''/'^-'/')du^ = -mlRl,u'/'^-y% (24) 
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or in its canonical form 

dl + 



,12 15 13 



duVn = -^^^^en (25) 



For Uq » Ma, /(m) — > 1, the e.o.m simplifies and the qualitative behavior of m„ can 
be determined by using dimensional arguments [13] . Define the dimensionless parameter 
V = ^ then for the limit Mq >> ma where / ^ 1 



The e.o.m in terms of v reads 



7 - ^ r (26) 



a.(.;^^/^-^/^)a„a = -ml^v'/'j-'/% (27) 



Since the L.H.S is dimensionless so must be the R.H.S and hence 

,2 _ ^ 

Using the relation ([9]) between uq and L we find 



(28) 



m„ cx (29) 

while the mass of the glueball is related to m^f, ~ For the case ua = uq^ L = irR so 
the glueball and mesons masses have the same scale. However in the general case where 
Uq > ua there are two different scales m„ ~ x ^ ~ "^s^- 

In order to find the exact spectrum of the eigenvalues of (125|) one can use the shoot- 
ing technic which is implemented by demanding symmetric or anti-symmetric boundary 
condition to ^„ at m = Uq and integrating the equation up to the u » Uq region where 
the solution could be matched to its normalizable asymptotic expansion. Of course this 
matching is only possible when the correct eigenvalues are being used and so one shoots 
with different eigenvalues until a matching is obtained . However there is a problem with 
these coordinates at n = uq since, -^^\u=uo oo ( see eq. ([6])). An odd perturbation 
to the classical configuration will cause no change in the shape of this singularity but an 
even one will, and so will also have a singular derivative. 

This problem is reflected in the singularity of the e.o.m (l25l) at u ^ uq. To see this 
behavior explicitly we change coordinate to a dimensionless parameter z as follows 

= ul + u\z^ (30) 

the eigenvalue problem (!25|) then becomes 



q2 ^ 5<2 1 iu\z 



+ uXz"^ z {u\\ u\z'^YI'^^\ " (Mq -I- m|2;2)7/3 



a45„=-;5?^^«„ (31) 
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where 7' stands for the derivative of 7 with respect to u. Since 

_ ul{8uf5ul) _ , _ uli8uf5ul) 
Iz^O - , Iz^O l-^^j 

we find that this equation has a regular singularity at 2; = 0! 

Indeed it was already noticed in [6] that by employing the 'shooting' technic only half 
of the spectrum could be found, namely only the odd modes where seen while the even 
ones could not be obtained, these modes that should had been obtained by setting the 
boundary conditions to 

^4z = 0) = l ; dMz = 0) = 0. (33) 

turned to be singular and could not be integrated. In [1] only the special case of uq = u\ 
was analyzed, in this case since lim^^^^^ duXd = 0, a smooth and nonsingular transfor- 
mation into cartesian coordinates is allowed via 

= u\ + u\{z^ + y^) ; X4 = i?arctan(-) (34) 
The corresponding action for y is (after setting u\ = 1) 

5 ~ / d'^xdz + u(zy(d,y? + 2yA (35) 

J I u[z) J 

inserting the expantion y = J2n=i Vn{,x^)yn{z) the e.o.m for y^ is 

9.yn + Y^^d^yn - = (TT^^- (^6) 

which is non-singular. For the more general case of uq > we were not able to find a 
similar coordinate transformation and hence we follow a different approach desribed in 
the next section. 



3 A regular e.o.m for the scalar fluctuation at the 
low temperature phase 

As we have seen above, we could not obtain the even modes of the fluctuation |^ around 
the classical curve because -^f- diverges at u = Uq. The issue of choosing a direction 
along which one should anlayze the fluctuations, has been discussed in the context of the 
stringy description of the Wilson like [21] . It was found that the safest approach is to use 

^If the classical curve X4^ci was odd, then the odd mode would become singular. 
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the fluctuation in the direction which is normal to the classical configuration. For our 
case the normal to the classical configuration at the tip u = Uq is along the u direction. 
Thus from here on we study the fiuctuation in the u direction, that is 

U(X4, X'') = Ucl{Xi) + ^(X4, x'') (37) 

our classical curve would be Uci{x4^) and as can be seen from (ED we have ^|^=o = 
so the point u{x4 = 0) = Uq pause no problem now! The quadratic action for these 
fiuctuation is (after setting u\ = 1) 



(lln^^ + 18ao + 3nii - I2u^ - 27ao{u^ + u^) - 2u^] 



^1 



2u^^P 

where Oq = u^f{uQ) and it should be understood that u = Uci{x4) and its formal expres- 
sion is 



u{x,)= dx,f{u){-—f'\ 



after plugging a mode expansion the e.o.m in its canonical form is 



(39) 



d'jn - {- + ^M^n-^^^^n (40) 

u uj ao 



[llu^^ + 18ao + 6nii - 12u^ - 27aoiu^ + u^) - 2u^) 



2aoU^ 

where Ux = dx^Ud- Since there is no analytic expression for the integral in ( l39ll . we 
obtained u{xi) numerically during the integration of eq. (HOj) when 'shooting' to find the 
eigenvalues of (HOj) . 

The resulted spectra are summarized in figures ([T]), ([2]) and (j3]). The following properties 
charcterize these spectra 

• The first observation one can make is that for Uq = u/^ our results for the symmetric 
and anti-symmetric lowest lying states match those of [1]. 

ml = 3.3 ■ ml = 5.3 (41) 

• The figures ([1]) and ([2]) describe the dependence of the squared mass of the first 
excited symmetric and anti-symmetric states as a function of the "constituent quark 
mass" defined in [5] and as follows 



1 fuo T /-MO 

~ 2^' Jus " 2'Ka' 



m': = : / ^-guQuudu = - — - / / ^''^{u)du (42) 



MA 
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This parameter relates to the constitutent quark mass and not to the current algebra 
( QCD) mass, since even when it is turned on the fluctuations that correspond to the 
pions are massless. In fact the quantity dual of the constituent quark mass should 
associate with plus a constant term which is independent of uq since already 
for Uq = Ma the mesons are massive and hence there is a non-trivial constituent 
mass. This assignment is also in accordance with the semi-classical description of 
high spin mesons [1] and their stringy split into two lower mass mesons [22] . /^From 
these figures we see that indeed for for uq > u\ the square of the mass of the scalars 
grows linearly with m^. This is to be contrasted with the results found in [5] for 
vector mesons of non-critical models where the masss itself is found to be linear 
with the TTig ( see also down in section (5).) 

• We have also determined the spectrum of the higher excited mesons both the sym- 
metric as well as the anti-symmetric ones. The dependence of the squared masses 
on the excitation number for various values of is drawn in figure ([3]). The linear 
fit to this curves are given by 

ml = 3.3 + 1.6r2^-^ ml = (43) 

ml = 10.5 + 6.5r2^-^^^ = 9.3 

ml = 15.8 + 9.5n^-^^^ m^ = 14.3 

Stringy modes are characterized by the well known ~ n behavior. We thus see 
that the scalar meson spectra that follow from the model of [T] do not correspond 
to stringy modes. This is of course of no surprise since it follows from a low energy 
effective field theory and not from a semi-classical treatment. 

• Last by not least we see from figure ([1]) that the lowest scalar excitation remain 
non-tachyonic for all values of Mq which serves as an partial evidence for the stability 
of the Sakai Sugimoto model. 



4 Scalar mesons in a non critical holographic model 

We would like now to find the masses of the scalar modes associated with the fluctuations 
of the probe brane around the classical profile in the non-critical gravity background of 
[7j. Using the background (fTS!) in an effective action that includes only the DBI. The CS 
term 

Scs-^ f C,^ [ (44) 
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The first symmetric mode ffi^vs m 
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The first antisymmetric mode irP vsm 
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(B) 



Figure 1: (A) The mass squared ml of the lowest exited symmetric mode as a function 
of {Roi = ua = 1) 

Figure 2: (B) The mass squared m2 of the lowest exited antisymmetric mode as a function 
of {Rda = ua = 1) 

Mass squared f the mesons vs. their excitationnumber at T<7'j for 
different values of m 
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Figure 3: (A) The tower of the mesons squared mass in the low temperature phase 
(i?z)4 = ut = 1) 
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would contribute upon a substitution of the form (1371) for u and expanding a linear and 
quadratic term to the action, thus affecting the spectrum. However, it was found in [8] 
that including this CS does not yield a sensible thermal phase diagram and hence we 
discuss separately an effective action that includes only a DBI action and one with both 
the DBI and CS terms. We start first with the former case: Analyzing the spectrum in a 
similar manner to the analysis of section (3) we find that the fluctuations are subjected 
to the following eigenvalue equation 

9„(A-3/^)a.en = -^^^^0^^n (45) 

Like in the critical case, for uq » the qualitative behavior of m'^ can be seen by 
changing the variable u into the dimensionless parameter y = —. At the limit Uo >> 
we find that f{u)^l and so 



and find that in terms of the dimensionless parameter y the e.o.m is now 

dy{y'r"')dyin = "^^^Cn (47) 

"0 7 

Since the L.H.S is dimensionless so is the R.H.S and we find 

< oc ^ (48) 

Note that due to the different background now L ~ and hence again we get that 
m'^ ~ j^. However in terms of the asymptotic behavior is that m'^ ~ m'^ and not 
m'^ ~ as was the case for the mesons of the critical model. 

Repeating the exact same steps as for the critical case we find that the quadratic 
action for fluctuation in the direction around the classical curve leads to an e.o.m 
which is singular at u = uq and as a consequence the attempt carried in [S] to obtain the 
spectrum of the even modes had indeed failed. And so like in the critical case we turn to 
study the fluctuation in the u direction instead. The action for the fluctuation is then 

1 r 1/2 p4 

al^\u^ + 36ao - 63m^° + Uu^° + 48m^^ - 9200^^ - 44aJ°) 

2m22/3 

and indeed this action leads to a regular e.o.m at m(x4 = 0) = uq. 

14 1^ 7/10 f 2 d4 

dl^ - (- + + ^^-^-^v'-d.d.^ 

u u^j ao 
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{u^ + 36ao - 63^1° + Mm^o + 43^15 _ 92aou^ - 44aoM^°) 
+ 2^0 ^" ^ ^ 

Using the shooting technic we found the eigenvalues of different modes of the fluctuation 
for various values of m^, our finding are summarize in figures (jl]),©. One can see that 
the masses m'l and m'g grow linearly with as expected from ()18|1 . At Uq = u\ = 1 we 
fincfl. 

m'/ = 1.51 ; mZ = 2.07. (51) 

which is in agrement with [5]. Again we also studied the dependence of the mass on 
the excitation number, the results are summarized in figure (jS]) and are: 

m„ = 1.51 + 2.32n^-°^ = (52) 

rUn = 13.5 + 4.95n^-°^ = 9.3 

Thus we see that both in terms of the dependence on n as well as the dependence on 
the scalar meson spectra admits a different behavior than that of the critical model 
of [T]. A similar behavior has been observed for the vector mesons in 
Next we consider the case where the effective action includes both the DBI and CS 
terms. Including now the CS term (with its full strengh a = 1) the quadratic action for 
the fluctuation becomes 

^ = I f HS^^-^"^^' ^ ^^i^^^"^^' 

51/2(^5 ^36B-63ni° + 14^20 + 48^15 -925^5 -4451°) 20 

— 11 f 

2n22/3 75 ^ 

where B = (mq/^/^(mo) — Mq + m^)^ and the e.o.m is then 

(u^ + 36B - 63ni° + Wn^o + 48u^^ - 92Bu^ - UBu^^) 20u^ , , 

+ ^+v!537^« = °(5^' 

With the Chers Simon taking into account the dependence of the mass squared on 
the excitation number is now to be read from figure ([9]) to be: 

m„ = 2.07 + 5.42n°-^^ = (55) 

m„ = 16.49 + 1.01ra°-^^ m\ = 9.3 

The dependance on m'l is described in figures ([7]) and ([8]). 



'^Our results are for RAdS = 1- 

®To keep contact with the resuhs in [5] we had rcnormahzcd the masses by the factor | coming from 
the change of variables u — > z. 
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Non critical model: Mass of the symmetric mode : m' vs tn 



Non-critical model: Mass of the antisymmetric mode: ni vs m 
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Figure 4: (A) The mass m'l of the lowest exited symmetric mode of the non-critical 
model as a function of {RAds = = !)• 

Figure 5: (B) The mass m'2 of the lowest exited antisymmetric mode of the non-critical 
model as a function of {Racis = = !)• 
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Figure 6: The tower of mesons masses in the non-critical model 
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Non critical model with CS term (a-l): The mass of the first excited 
symmetric mode vs. m . 



M, 



Non critical modelwith CS term (a=l) : The mass of the first 
excited anti-synmietric mode vs. tn . 
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Figure 7: (A) The mass m'l of the lowest exited symmetric mode of the non-critical 



model with CS term included function of m'i . 



Figure 8: (B) The mass m'2 of the lowest exited antisymmetric mode of the non-critical 



model with CS term included clS cl function of -nf;, . 



The non critical model with CS term (a-l): The tower of the mesons 
mass vs. their excitation number. 
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Figure 9: The tower of mesons masses m'„ in the non-critical model with CS term included 



16 



5 Parity and charge conjucation 



In order to compare the resulting spectra from both the critical and non-critical models, 
we first have to identify the "quantum numbers" of the states that correspond to the 
fluctuations. More explicitly we have to determine the operations in the gravity mod- 
els which correspond to charge conjugation and parity transformations. In the model 
of [l] they were defined as follows: The charge conjugation operation associates with 
exchanging the left and right handed quarks which maps into the interchange of a D8 
and an anti D8 or differently transforming z — > —z. Parity transformation in the five- 
dimensional space-time spanned by Xi,z where i = 1,2,3 means the following transfor- 
mation {xi,z) {—Xi,—z). 

For the generalized set up with uq > u\ we can still define the coordinate z as follows 

= ul + uaz"^ (56) 

Note the difference with respect to ( l30l) since here we take z to have dimension of length. 
With this definition of the z coordinate the discrete transformations of [T] remain in tact. 
The effective action on the probe brane has to be invariant under both parity and charge 
conjugation. The DBI part (12T!) is quadratic in ^ and hence cannot determine the right 
transformation of the fluctuation modes. The situation with the CS term is different. 
Recall that the CS term has the form 

5c5~/ FAFAC5=f F AFA [ (fxdzC^ = ( F AF ( d^xdz^{x^, z) (57) 

The last part we have used the explicit form of the C5 

C5 = ^(x'^, z)dx^ A ...dx^ A dz (58) 

In order for this term in the action to be invariant under parity and charge conjugation 
it is clear that ^{x^, z) has to be even under both charge confugation and parity trans- 
formation. Now since ^(x, z) = Y.n fn{,x^)in{,z) we conclude that the map between the 
fluctuation modes and sclar particles is the following 

symmetric C,n O"*"^ mesons 

antisymmetic C,n mesons (59) 

For the non-critical model again the DBI action does not determine the transforma- 
tions of ^ under parity and charge conjugation. We have argued above based on [8] that 
a CS term of the form (14^ should not be incorporated. Thus there is no way to this 
order to determine the transformation of ^. 
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Without the constraint from the CS term we may have that ^ is even or odd under 
charge conjugation and parity transformtaions. In the latter case the assignments of (!59ll 
have to be reversed, namely symmetric ^ corresponds to and antisymmetric ^ to 0''""^ 

Next we want to compare the spectra to mesons observed in nature. It is well know 
that scalar mesons in nature are either 0"'"'' or pseudo sclars of the form ^ and there are 
no observed low lying mesons of the form . Thus there is a serious mismatch between 
the holographic scalar mesons extracted from models with flavor branes anti-branes of 
critical models and with the observed mesons in nature. We will come back to this issue 
in the conclusions. 



6 Scalar mesons in the intermediate temperature phase 

The background that corresponds to the deconfined phase namely T > 1/2ttR is given in 
(fTT]) . As was shown in [13] this deconfined background can admit also a phase where chiral 
symmetry is broken, the so called "intemediate phase" We now analyze the spectrum of 
the scalar mesons in this phase. Since the procedure of extracting the scalar meson is 
identical to that of the low temperature analysis of the previous sections we present the 
final results for the spectra of masses. The spectra are presented in figures ( fTOl) . ( fTTj) and 
([T2|) . The main features that these spectra admit are the following 

• As can be seen, at the phase transition T = Td the value of the masses are (for the 
values My = 1, Mo = 8) 

m^(T = Trf) = 8.36 ; mL(T = T^) = 45.96 (60) 

while in the low temperature phase at the point of phase transition with u\ — >■ 
ut = 1, uq = 8 the masses are 

m^(T = Trf) = 8.40 ; m2^(T = T,) = 46.00 (61) 

We see a very small jump in the masses at the transition point, the same as was 
seen for the vectors in [6] 

• While in the low temperature confining phase the masses of the mesons are tem- 
perature independent since the background in this phase does not depend on the 
temperature, the masses of the mesons do depend on the temperture in the inter- 
mediate deconfined phase. As was observed in Lattice simulations and was found 
also for holographic vector mesons [6], the masses decrease as a function of the tem- 
perature. The symmetric mesons decrease at the chiral symmetry phase transition 
temperature T = T^sb to ~ 60% persent of their values whereas the antisymmetric 
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ones to ~ 80%. This drop off is much more significant than for the vecotr mesons 
of the critical model [6]. 

• Note that it is only consistent to increase the temperature up to where the next 
phase transition occur and chiral symmetry is restored. 

This happens at T = T^sb (for the choice mq = 8 we found that T^sb = 2.44Trf), 
then the merged together D8-D8 breaks into a separate pair of D8-D8. We can 
also see from figure ( fTOl) that if we continue to increase the temperature beyond 
T^SB then at some point the scalar mode becomes Tachyonic, signaling that this 
background is no longer stable at this temperature as indeed we know. 



Like in the low temperature we had also checked the squared masses dependence 
on the excitation number see figure (fT2|) . this was found to be: 



ml = 8.3 + QAn^-^ T = (62) 
ml = 7.6 + 6.9n^-^^ T = 2Td 



7 Conclusions 

I this paper we had overcome technical problems faced in [221 E] and succeeded to obtain 
the holographic mass spectra of the scalars in the low and intermediate phases of the 
chiral symmetry broken phase of the critical model and also of those of the non-critical. 
Let us summarize the results of this work and mention certain open directions. 

• There is a difference between the dependence of the mass of the scalar mesons on 
the "constituent mass parameter" m^. In the ten dimensional models one finds a 

oc relation (see figures ([I]),® for the first two excited modes), whereas for 
the non-critical model the relation is m oc (see figures (I1D,([S]) and ([T]),®)- 

• Both the critical models and the non-critical one do not admit a Regge/stringy 
behavior of ~ n. This is not unexpected since the stringy excitations is not 
visible in the low energy effective field thoery. 

• One can compare the ratio of the low lying mesons both vector and scalar mesons 
to those observed in nature. Table (1) present such a comparison. It is interesting 
to note that turning on a constituent mass improves the ratios with respect to 
those for zero m'i. 
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Figure 10: (A) The mass squared m\{T) of the lowest exited symmetric mode as a 
function of T/T^ (mq = 8, Roi = 1 and R = 2/3) 

Figure 11: (B) The mass squared niKT) of the lowest exited antisymmetric mode as a 
function of T/T^ {uq = 8, Rm = 1 and R = 2/3) 
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Figure 12: The tower of mesons squared mass in the intermediate phase 
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Table 1: A comparison to the experimental data where the best fitted is presented 
vs. mjj = (for the critical case we have found that there is no improvement in ratios of 
the vectors so we have left these entries empty.). 





experiment 


D4-D8 at = / 0.38 


Non-critical at = / 0.16 




2.51 


2.4 / - 


2.8 / 2.62 




3.56 


4.3 / - 


5.5 / 5.29 




3.61 


4.9 / 3.63 


4.1 / 3.65 




0.7 


0.49 / 0.62 


0.67 / 0.75 



• The hologrphic spectra of the critical models admit a branch of scalar mesons of the 
type . These does not exist in nature. It seems to be a severe shortcoming of 
these holographic models. This difference cannot be attributed to the fact that we 
consider large iV^. It will be interesting to investigate the question of how generic 
is this situation and whether one can construct a mechanism to project it out from 
the low lying spectra. 

• The behavior of the scalar mesons at finite temperature in the intermediate phase 
is similar to that for the vector meson in the model of [6] . However the decrease of 
the mass with increasing temperature is more dramatic for the scalar mesons. It is 
interesting to check if a similar phenomenon occurs also in lattice simulations. 
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